Competition between Multiple Totally Asymmetric Simple Exclusion Processes for a 

Finite Pool of Resources 
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Using Monte Carlo simulations and a domain wall theory, we discuss the effect of couphng several 
totally asymmetric simple exclusion processes (TASEPs) to a finite reservoir of particles. This 
simple model mimics directed biological transport processes in the presence of finite resources, such 
as protein synthesis limited by a finite pool of ribosomes. If all TASEPs have equal length, we find 
behavior which is analogous to a single TASEP coupled to a finite pool. For the more generic case 
of chains with different lengths, several unanticipated new regimes emerge. A generalized domain 
wall theory captures our findings in good agreement with simulation results. 



I. INTRODUCTION 



A fundamental and comprehensive understanding of 
non-equilibrium phenomena remains one of the great- 
est challenges of current condensed matter and mate- 
rials physics [H, with significant consequences for ad- 
vances in materials science, the life sciences, and en- 
gineering. Even the non-equilibrium steady states of 
open, transport-carrying systems continue to defy our 
equilibrium-trained expectations and intuitions. One ap- 
proach towards progress has focused on investigating sim- 
ple model systems, with the goal of identifying generic 
classes of behaviors. 

The totally asymmetric simple exclusion process 
(TASEP) 0, [1| is one of these models. It has ac- 
quired paradigmatic status for several reasons: (i) it is 
very simple, consisting of particles hopping along a one- 
dimensional chain; (ii) with open boundaries, it shows 
highly nontrivial behaviors, such as distinct phases, 
shocks, and long-range correlations in both space and 
time; (iii) in its simplest forms, its steady-state proper- 
ties, as well as selected dynamic quantities, can be found 
exactly; and finally, (iv) the model is closely related to 
interesting applications, such as biological transport [3| 
or traffic flow [5]. At the origin of this richness lies the 
violation of detailed balance. The specific behaviors de- 
pend strongly on the boundary conditions. With peri- 
odic boundary conditions, the stationary state is a flat 
distribution with all configurations equally probable Q. 
However, the dynamics of this system is nontrivial and 
differs from simple diffusion Q. With open boundary 
conditions, particles are injected at one end and removed 
at the other with different (but constant) rates. In this 
case, even the steady state is nontrivial and remained 
unknown for two decades 0, [l^, [llj. Despite being a 
one-dimensional system with short range interactions and 
dynamics, the open TASEP displays distinct (stationary) 
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phases Q, controlled by the entrance and exit rates. As 
may be expected, the dynamic properties are even more 
complex and rich [H [ll, [ll [l3 . 

These TASEP studies have recently been extended by 
coupling the chain to a finite (rather than infinite) parti- 
cle reservoir [l^ , reflecting a constraint on the total 
number of particles. In [l^, the particles represent cars 
leaving a parking garage, so that the rate of entry onto 
the roadway (the lattice) is chosen to be a constant, as 
long as there is at least one car in the garage. In [17l |. 
the TASEP models a biological transport process 0] , and 
the constraint reflects the finite number of ribosomes in 
a cell, with those leaving the mRNA (the lattice) being 
"recycled" to the entry point. The origins and effects of 
"ribosome recycling" are multifaceted, such as diffusion 
of the recycled components from termination to initiation 
sites [l^. Addressing all possible issues for a real cell is 
beyond the scope of this study and our aim here is mod- 
est, namely, to explore how the finite pool of available 
particles affects the lattice occupation and currents. As 
in , we consider an entry rate that depends smoothly 
on the number of particles in the pool, Np. In particu- 
lar, we choose a rate function which is proportional to 
Np when the concentrations of ribosomes in the cell are 
limited, and, when particles are abundant, becomes a 
constant - the inherent binding rate of a ribosome. The 
simulation results for a single TASEP can be described 
well by a domain wall (DW) theory [isl. [l^ [l7| , esp ecially 
when generalized to account for the feedback [I^l . 

In addition to being constrained to a finite pool of ri- 
bosomes in a living cell, a mRNA must compete against 
many other different genes (or mRNAs) for these re- 
sources. Therefore, we are motivated to study the com- 
petition of multiple TASEPs for the same pool of parti- 
cles. Since different mRNAs are comprised of different 
numbers of codons, it is reasonable to study TASEPs 
with different lengths. Does one "win" and another one 
"lose"? What does "winning" and "losing" mean? 

This paper is organized as follows. The next section 
introduces our model. In Section IIIIl we present simu- 
lation results for two and three TASEPs connected to a 
finite reservoir of particles. Analytic results, based on our 
generalized DW theory [l^, are derived for an arbitrary 
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collection of TASEPs and will be discussed in Section lTVl 
A final section is devoted to a summary and outlook. 



II. A MODEL FOR COMPETING TASEPS 

Let us first review the standard (or "ordinary") 
TASEP. Each site of a discrete lattice of length L, is 
labeled by i = 1,...,L and may be empty or occupied 
by a single particle. Thus, a particular configuration can 
be described by a set of occupation numbers, {n{i)}, 
each taking the values 0, 1. A configuration evolves to 
a new one according to the following rules. A particle 
is chosen at random, say, at site i. If the neighboring 
site to its right (site i + 1) is empty, the particle hops 
there with rate unity. If the particle is located on the 
last site (i = L) of the lattice, it exits with a proba- 
bility p. In addition to the particles on the lattice, we 
assign a "virtual particle" to label an external reservoir, 
so that when chosen, a particle is placed with probabil- 
ity a, on the first site {i = 1) of the lattice, provided 
n (1) = 0. Notice that this system can be regarded as a 
total of Ntot particles hopping on a ring with L+1 sites, 
where the site i = is associated with the reservoir and 
can be occupied by any number of particles. Of course, 
the hopping rates into and out of this site are different 
from the rest. They can be represented, respectively, as 
l3n [L) and a[l — n {!)]. Note that neither depends on 
n (0) and so, as long as Ntot > L , there would be at least 
one particle which can be injected onto the lattice. In- 
deed, this is the scenario presented in the parking garage 
problem [l6| , in which interesting transitions are studied 
for Ntot G [0,L]. Lastly, the seemingly complicated rule 
of choosing real and virtual particles can be replaced, in 
this formulation of TASEP, by: "Randomly choose an 
occupied site." 

In the steady state, this simple model exhibits three 
phases: a low density (LD) phase for a < 1/2 and a < 
(3, a high density (HD) phase for (3 < 1/2 and (3 < a, 
and a maximal current (MC) phase for a, (3 > 1/2. The 
densities p in the three phases are given by a (LD), 1 — P 
(HD), and 1/2 (MC), respectively The line a = /3 < 
1/2 marks the coexistence of HD and LD phases and is 
characterized by the presence of a shock which performs 
an unbiased random walk through the whole lattice. For 
this reason, the coexistence line is also sometimes referred 
to as the "shock phase" (SP). Since the entrance and exit 
rates are constant here, Ntot plays no role and can be 
finite (but larger than L) instead of oo. 

Our motivation for studying the TASEP is modeling 
protein synthesis in a cell [J], in which the particles rep- 
resent ribosomes. Thus, Ntot is finite and must be shared 
between many mRNA's (there may be as many as 10,000 
mRNAs in some cells). Only the unbound ribosomes (i.e., 
those in the "pool", totalling Np) are available for initi- 
ation (i.e., to enter any one TASEP). Assuming the con- 
centration of these ribosomes is uniform, we introduced 
a more realistic model flf\ for the entry rate, a, that 



depends on Np as follows. Starting with just one lattice 
(mRNA) in our system, let us denote its total occupation 
by A^, so that 

Ntot = N + Np (1) 

is kept fixed. Particles still exit the lattice as before, 
with rate (3, and are placed immediately into the pool. 
In contrast to the ordinary TASEP, the (effective) en- 
try rate, cteff, will now be controlled by Np, through 
aef f = af{Np). Here, a is a constant and f{Np) assumes 
a value in [0, 1] for all Np. Physically, a^ff should vanish 
if there are no particles in the pool, whence we impose 
/(O) — 0. Further, / should be a monotonically increas- 
ing smooth function of Np. Finally, for sufhcicntly large 
Np, we should recover the standard TASEP with a con- 
stant entry rate, a, whence liuiN^^oo f (Np) = 1. These 
properties are motivated by the notion that the rate for 
a ribosome to bind to the mRNA (a) should be limited 
by the ribosome availability, especially for low concen- 
trations and (b) should approach some intrinsic rate for 
the binding process, when the ribosomes are abundant. 
The specific form of f{Np) is not very important but we 
follow 17] and use 

fiNp) = tanh(A^p/A^*) (2) 

where N* is a suitably defined normalization, or 
crossover, scale. In [l7|, N* was chosen to be the av- 
erage number of particles for the standard TASEP with 
entry and exit rates a and f3. 

In the following, we will consider multiple open 
TASEPs of various lengths, Li {i = 1,2,...M), M be- 
ing the total number of chains in our system. Writing 
the occupation in each as Ni, we naturally define 

Pi = Ne/Li (3) 

as the overall densities on each chain and of course, 

M 

Ntot = J2^f + ^P (4) 

e=i 

as the generalization of Eqn. ([T]). Most of our simulation 
results are for M = 2, with a few for Af = 3. All M 
TASEPs draw their particles from the same reservoir, 
according to the same cteff = o:f{Np) with / given by 
Eqn. For multiple chains of different lengths, we 

choose N* to be 

M 

N* = M-^J2piLi (5) 

where pi is the average density for an ordinary TASEP 
of length Li, with entry and exit rates a and This 
choice of N* only serves to locate the specific value of 
Ntot at which the system crosses a phase boundary. It 
also allows us to define (somewhat arbitrarily) 

M 

ptot = Ntot/M-^ J2 ■ (6) 
1=1 
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For our simulation studies, we define a Monte Carlo 
Step (MCS) as making M + ^ attempts to choose 
a site to update. Note that, on the average, the pool 
is updated M times as often as a site on any lattice. 
This choice corresponds to updating all links (nearest 
neighbor pairs of sites) with equal probability, with the 
understanding that each TASEP is linked to the pool 
(independently of the others). In this study, the entry 
rates to all chains are chosen to be the same aeff- We 
typically start with all particles in the pool and wait 100k 
MCS for the system to reach steady state. Every 100 
MCS thereafter, we measure the site occupations, n{i), 
in each chain. Using up to 10k measurements as our 
"ensemble," we compute the density profile 



p{i) = {n (i)) 



(7) 



and the overall density p = / ^'^^ each chain. 

The average steady-state current, denoted by J, is ob- 
tained by monitoring the total number of particles that 
enter (and exit) a chain over the last IM MCS. Note that 
J is also just the product Pp (L). 



III. SIMULATION RESULTS 

In this section, we present simulation data for the 
model described above. In principle, any number of 
TASEPs can be connected to the pool. However, to gain 
some insight into competition, we begin with the simplest 
case, with only two TASEPs. Although some selected re- 
sults for the M = 3 case will be presented at the end, 
most of our data here is for the case of M — 2. As 
we will see, novel features already appear when just one 
more TASEP is added to the system. By contrast, we 
have not encountered, so far, any further surprising phe- 
nomena by considering M > 2. As in the earlier study 
[l3|, we explore four regimes here, with a and (3 cor- 
responding to the LD, HD, MC phases of the ordinary 
TASEP, as well as the coexistence line, SP. Our main in- 
terest is how an increasing Ntot affect the various profiles 
and so, the average overall densities and currents. 



A. LD phase 

When a < 1/2 and a < f3, the standard TASEP dis- 
plays the LD phase. If coupled to a finite pool, it remains 
in the LD phase, since aeff < a. The only difference 
is that the density and current of the system rise lin- 
early with Ntot for small values of Ntot, and approach 
their asymptotic values from below as Not — > oo. When 
two TASEPs compete for this finite pool of resources, 
similar behavior is found. The results for equal length 
chains are illustrated in Fig. |l(a) for a — 0.3, /3 ~ 0.7, 
and Li ^ L2 = 1000. As Ntot is increased, the two 
TASEPs both remain in the LD phase, with equal den- 
sities, pi = P2 = P, and currents, Ji = J2 = J ■ Since 



both TASEPs are controlled by the same aeff and /9, 
there is full symmetry between the two and the observed 
behavior is hardly remarkable. The only notable differ- 
ence between our system and the one with a single con- 
strained TASEP is observed in p{Ntot) for Ntot <C aL. 
For our case, p — > ptot/i as opposed to ptot/'i for the sin- 
gle TASEP. This property is easily understandable, since 
we can regard the pool as an additional "chain" and note 
that the resources are evenly distributed amongst three 
(or two) "chains." For large Ntot, both TASEP den- 
sities level off to a, of course. Slightly more interest- 
ing is the case where the two chains have very different 
lengths, such as Li = 1000, L2 = 100. Due to the lack 
of symmetry, it is less obvious that pi — P2 ^ Ptot/i 
is still valid for small ptot a. However, it is straight- 
forward to show, using the methods in previous studies 
[13, [Hi, that pi = ... = pm Ptot/ (M + 1) in general, 
given the specific form of aeff {Np) we chose. At the 
opposite limit, the approach to the asymptotic regime 
is somewhat faster than for the equal length case. This 
behavior is also expected, since the "total" system size 
(Li + L2 — 1100) is considerably smaller than the case 
above (Li + L2 — 2000), so that saturation sets in at 
smaller values of Ntot- 

To summarize, the overall densities of competing 
chains behave just as a single TASEP coupled to a finite 
pool of particles. Meanwhile, the overall currents are, 
except for finite size effects, also essentially the same: 
P(l-P)- 



B. MC phase 



For a single TASEP coupled to a finite pool with 
a,(3 > 1/2, the current approaches its asymptotic value 
(1/4) smoothly from below. However, the density ex- 
hibits a sharp "kink," which marks the crossing of the 
LD-MC phase boundary when Ntot becomes large enough 
to sustain a density of 1/2 on the chain. The range of 
Ntot over which this crossover occurs is controlled by 
the system size: For large system sizes, it becomes very 
narrow and therefore, appears as a "kink" ; in smaller 
systems, the crossover occurs over a larger range and 
appears smooth. For the case of two TASEPs, we ob- 
serve very similar behavior. Fig. [2] shows our data for 
a = P = 0.7 and Li = L2 = 1000. Again, the currents 
and densities of the two TASEPs coincide, provided the 
two TASEPs have equal lengths. For unequal lengths, we 
observe the anticipated finite-size effect in the density: 
The longer TASEP displays a much sharper crossover 
from LD to MC behavior than the shorter one. Like the 
LD case, the overall densities and currents of two com- 
peting chains behave much like those for a single TASEP, 
including the finite size effects. 
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FIG. 1: (Color online) Average overall density and current as 
a function of Ntot for [(a)] Li = L2 = 1000 and [(b)] Li = 1000 
and L2 = 100 with a = 0.3 and P = 0.7. 
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FIG. 2: (Color online) Average overall density and current as 
a function of Ntot for [(a)] Li =1/2 = 1000 and [(b)] Li = 1000 



and 1/2 = 100 with a = /3 = 0.7, 



C. SP 

The SP case is the most challenging, due to the large 
fluctuations that characterize this "phase" . In the stan- 
dard TASEP, the system exhibits a freely moving shock, 
separating a low-density from a high-density region. In 
the constrained TASEP, the crossover from the LD to the 
SP phase is highly nontrivial. The details depend on both 
the length of the chain and the functional form of f . 
Fortunately, the essentials are well captured by DW the- 
ory |13|.[l4l.ll7l|. especially when it is generalized appropri- 
ately [1^1 . Not surprisingly, two TASEPs of equal length 
exhibit the same densities and currents as a function of 
Ntot- Differences emerge only when Li ^ L2, as illus- 
trated in Fig. m where a = /3 = 0.3 and Li = 1000, 
L2 = 100. These differences may be expected, however, 
if we recall that single TASEPs with various lengths be- 
have quite differently when coupled to a finite pool [13] ■ 
Details of the origins of these differences in our case of 
two competing TASEPs can be understood in terms of 
the generalized DW theory fl9!| which we discuss in the 



next Section. Again, like the cases above, when finite size 
effects are taken into account, there are no qualitatively 
new phenomena when an additional chains is included in 
the competition for a finite pool of particles. 



D. HD phase 

The most interesting results are observed with a and 
/3 set in the HD phase. For a single constrained TASEP, 
there are three distinct regimes in the average total den- 
sity as Ntot is varied. The behaviors for small and 
large Ntot are expected: For the former, they follow the 
LD-behavior discussed above and for the latter, the den- 
sity and current approach their asymptotic values. The 
presence of the third regime, interpolating between these 
two limits, was a surprise initially. Here, the density rises 
linearly with Ntot while the current remains constant at 
its asymptotic limit. The origin of this linear depen- 
dence can be traced to Np (the reservoir occupation) be- 
ing essentially constant over a range of Ntot , so that any 
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FIG. 3: (Color online) Average overall density and current as 
a function of Ntot for [(a)] Li = L2 = 1000 and [(b)] Li = 1000 
and L2 = 100 with a = /3 = 0.3. 



changes in Ntot are absorbed by the lattice. In particular, 
Np takes the critical value N^, given by ttg// (Np) — f3, 
so that coexistence of low and high density regions on 
the lattice is possible. Unlike the ordinary TASEP how- 
ever, the shock does not wander over the entire lattice. 
Instead, it is localized at some point determined by Ntot, 
with an intrinsic width controlled by d^^f [l^. As for 
the current, it has already reached /3(1 — /3) at the lower 
extreme of this linear regime and so, it is not sensitive to 
the final crossover into the asymptotic HD regime. 

Turning to a system with two TASEPs, there appear 
to be no new surprises if they are of equal lengths, as dis- 
played in Fig. [Hfor Li = L2 = 1000. The densities and 
currents for both TASEPs coincide and exhibit the three 
distinct regimes found in the single TASEP case. In stark 
contrast, however, a new feature emerges if Li ^ L2. To 
be specific, we will discuss the case Li = 1000, L2 — 100 
here, unless otherwise stated. While the density of the 
longer TASEP still displays the expected three regimes, 
the shorter system now develops five regimes, as illus- 
trated in Fig. [51 The two outermost regimes (for Ntot 
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FIG. 4: (Color online) Average overall density and current 
as a function of Ntot for Li = L2 = 1000 with a = 0.7 and 
/3 = 0.3. 
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FIG. 5: (Color online) Average overall density and current as 
a function of Ntot for Li = 1000 and L2 = 100 with a = 0.7 
and /3 = 0.3. 



near zero and Ntot — *■ 00) are the familiar LD and asymp- 
totic regime. In the central region, however, we see a 
horizontal plateau, bordered by two "crossover" regimes 
where the density increases steeply with Ntot- The cur- 
rents show no discernible signature whatsoever, exhibit- 
ing only the LD and the saturation regime. To appre- 
ciate the different behaviors in intuitive terms, we will 
present several perspectives here. In the next section, 
we will provide the mathematical approach, in which an 
exact solution of the generalized DW theory appears to 
account for all data quite well. 

First, let us examine the "plateau" region. Similar to 
the case of a single constrained TASEP, this region is 
characterized by Np being fixed at N^, and so, cteff — P 
while changes of Ntot are absorbed by the two TASEPs. 
The main difference here is that the excess particles are 
free to choose either chain. Given that the pool suffers 
only minor fluctuations, A^i -|- N2 is also relatively con- 
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stant, so that the two chains simply trade particles back 
and forth. In terms of the domain wall picture, a DW 
is present in each lattice, but their walks are completely 
anticorrelated. Though the shocks are no longer local- 
ized as before ^9] and free to wander about, they are 
limited by (Li,L2). Specifically, the shorter TASEP im- 
poses the extent by which the DW may wander on the 
longer lattice. In the "plateau" region, the DW is free 
to be anywhere on the shorter lattice, so that the av- 
erage profile there is strictly linear. The overall density 
here is strictly 1/2 and the "plateau" emerges. In this 
region, all changes of Ntot are absorbed entirely by the 
longer TASEP (i.e., dpi/ dN tot is precisely in Fig. 

[5|). However, as will be shown below, the similarity of 
this behavior to the single TASEP case is deceptive. The 
density profiles are quite distinct, reflecting the different 
origins of the constraints on the DW walks. This picture 
also shows clearly why such a "plateau" cannot exist for 
the Li = L2 = L case: Neither TASEP imposes a limit 
on the other. By symmetry, the only point when a DW 
on one lattice can wander over one entire lattice is also 
the point where its counterpart can wander over all of 
the other lattice. At this point, both profiles must be 
linear and both densities are 1/2, so that TVi -I- A''2 = L. 
Since the pool occupation must remain at N^, this single 
point corresponds necessarily to Ntot = Np + L. 

Another perspective on the existence of the "plateau" 
is offered in Fig. [6l which shows schematic views of the 
domains in the N1-N2 plane in which we are likely to 
find the Li ~ L2 — 1000 system (Fig. and the 

Li = 1000, L2 = 100 case (Fig. Eb). For very small Ntot, 
both TASEPs will be in the LD phase, represented by the 
circles (green online) near the origin. As Ntot increases 
beyond this regime, the system will find itself mostly in 
the long ellipses (red online), aligned along lines of con- 
stant iVi -I- 7V2. Finally, for large Ntot, both lattices will 
saturate in the HD phase, indicated by the circle furthest 
from the origin (blue online). These two simple sketches 
bring out the main difference between the two cases: If 
the lattices are of unequal length, there is a range of Ntot 
for which the size of the ellipse remains the same. In 
this range, the average density on the shorter TASEP is 
about 1/2 (i.e., N2 = 100/2 here) while the occupation 
on the longer lattice continues to increase. When the 
data in Fig. [5] are replotted in the N1-N2 plane, this 
picture is largely confirmed: Fig. \7\ In contrast, for a 
system with Li — L2, there is no such comparable range 
and so, there is no "plateau" region. In Fig. [51 we pro- 
vide two examples of Monte Carlo data associated with 
these sketches. They are histograms for finding A^i and 
A'2 particles in the two TASEPs. The "ridges" in these 
plots correspond to two of the ellipses in Fig. [S] (a) along 
iVi iV2 = 1100 for the Li = L2 = 1000 case and (b) 
along N1 + N2 = 525 for the other case. These ridges also 
indicate that, subjected to N1+N2 = const., all (TVi, A^2) 
pairs are equally likely. Of course, these observations re- 
fiect just the picture presented above: the DWs on each 
chain are free to wander, but strongly anticorrelated and 
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FIG. 6: (Color online) Rough sketches of distributions in the 
plane for [(a)] Li = L2 = 1000 and [(b)] Li = 1000 and L2 = 
100. The dashed boxes correspond to the region allowed by 
LD/HD co-existence (i.e., densities of 0.3/0.7 here). 
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FIG. 7: (Color online) Simulation data for L\ — 1000, and 
L2 = 100 with a = 0.7 and /? = 0.3 in the N1-N2 plane. The 
dashed boxes correspond to the region allowed by LD/HD 
co-existence. 



limited by (Li, ^2)- 

To provide yet another perspective, we present the den- 
sity profiles of both chains in the Li = 1000, L2 = 100 
case in Fig. [9l Note that, in this figure, the horizontal 
axis is the fractional distance (i/L) along a chain and 
so, points on the two curves correspond to different ab- 
solute distances (i). Being linear, the profile of the short 
TASEP clearly reflects a totally delocahzed DW. For the 
long chain, though the profile indicates a localized shock, 
this appearance is deceptive. The width (over which 
the profile changes from LD to HD) is actually some- 
what larger than its counterpart in a single constrained 
TASEP (see, e.g.. Fig. 2 in Ref. The origin of 
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FIG. 8: (Color online) Distribution of iVi and N2 for [(aj 
Ntot = 1400 and Li = La = 1000 and [(b)] iVtot = 700, U = 
1000, and L2 = 100 with a = 0.7, and (3 = 0.3. The z-axis 
scale for P{Ni,N2) is arbitrary in both figures. 



this broadening can be traced to the additional fluctu- 
ations allowed by the shorter chain (approximately 100 
lattice sites in this case). In particular, we re-examined a 
single TASEP of 1000 sites at the same a and f3 and ob- 
tained its profile. To accentuate the "interface" of the 
profile, we plot in Fig. [TUTa) the local slopes of the 
profiles (i.e., Ap{i) = p{i + l) — p{i), shown as open 
diamonds). Note that these can be regarded roughly 
as the probability of finding the shock at site i. Next, 
we construct a "smeared" data set Ap{i\m) as follows: 
Shift the raw data by j = l_^.,m sites and then aver- 
age over these m data sets [21|. By this procedure, we 
hope to account for the extra degree of freedom which the 
shock experiences, thanks to the presence of the short 
TASEP. The lines (solid red, thin dashed green, short 
dashed blue) in Fig. [TUT a) illustrate the result of smear- 
ing with m = 100, 200, and 300. Now we return to the 
two TASEP system and, in Fig. [10i;b,c), display Ap{i) 
(solid squares) for the Li,L2 = 1000,100 and 1000,300 
cases, respectively. Meanwhile, the lines here (solid red 
and short dashed blue, respectively) are precisely those 
from the smeared profiles in Fig. [TOT a). The good quan- 
titative agreement between the data and these Ap{i\m) 
is a convincing confirmation of the picture we presented: 
By freely exchanging particles between the two chains. 
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FIG. 9: (Color online) Average density profile for Ntot = 700, 
Li = 1000, and L2 = 100 with a = 0.7 and /3 = 0.3. 



the shorter TASEP provides extra room for an other- 
wise "localized" DW to wander in the longer chain. Our 
conclusion here is that there are two contributors to the 
localization of the DW in a TASEP competing for finite 
resources. They are (i) the feedback [19] due to a non- 
trivial aef / , producing an "intrinsic" localization length, 
and (ii) the constraint from the other TASEP participat- 
ing in the give-and-take of particles. Clearly, (ii) means 
that the longer chain imposes no constraint on the shorter 
one, so that the shock is completely delocalized, regard- 
less of mechanism (i). By the same token, the profile of 
the longer TASEP is typically somewhat more complex, 
as both mechanisms play a role. Naturally, if the chains 
are of equal length, delocalization prevails on both and 
the profiles will be linear (to the extent allowed by Ntot), 
as our simulations confirm. Now, despite the effects of 
competition, it is possible to isolate the role of mech- 
anism (i) and observe an "intrinsic" profile, as follows. 
For each measurement of {ni (z)}, we use the totals iVi.2 
to estimate the position of the DWs in each lattice: ki. 
Then, we average the shifted occupations {rii (i — ki)} to 
arrive at the "intrinsic" profile. The result is statistically 
identical to the profile in the single TASEP case [l^, so 
that we are confident of the merits of the intuitive pic- 
ture presented here. Details of this procedure and the 
comparisons will be provided elsewhere [25| . 

Having addressed the central "plateau" regime, we 
turn to the two bordering, "crossover" regions. In Fig. [51 
we see that the behaviors displayed here are quite rich. 
To highlight these better, we plot the gradients of the 
three occupations, iVi.2,p7 associated with increasing Ntot 
by 10 particles. In Fig. [Til we show the data from the 
more interesting Li — 1000, L2 — 100 case. The re- 
gions corresponding to the central plateau and its bor- 
ders are most clearly seen for the small chain (A^2 , solid 
diamonds, red online): two peaks with a flat valley in 
between. At a very naive level, these features can be 
roughly understood from the sketch in Fig. |6(b)[ For 
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while A''! continues to increase. A similar crossover re- 
gion is present for the right end of the rectangle, when 
again increases. Though this intuitive reasoning provides 
a qualitatively picture of the five regimes, it clearly fails 
to capture the details of the two crossover regions. For 
example, the changes in the two -/V's in the first crossover 
are anticorrelated, rather than increasing together. Evi- 
dently, these details are sufficiently subtle that they can 
only be fully understood in a quantitative theory for do- 
main wall motion - the subject of Section HVl 



E. Competition between three TASEPs 
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FIG. 10: (Color online) (a) Local slopes of the profile for 
a single TASEP with L = 1000 and its "smeared" versions, 
shown here as lines (solid red, thin dashed green, short dashed 
blue). Open diamonds (0) are data points, obtained with 
a = 0.7, /? = 0.3, iV* = 385, iVtot = 650, and shifted so that 
the peak lies at the origin. (b),(c) Solid squares (■) are data 
points for the L — 1000 TASEP competing with a shorter 
one: L = 100 (b) and 300 (c), respectively. Lines are same as 
those in (a). 
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FIG. 11: (Color online) Occupation gradients of A'^i, N2, and 
Np with respect to Ntot for Li = 1000, and L2 = 100 with 
a = 0.7 and /3 = 0.3. 



small Ntot, the circles (green online) traverse along a 
line with slope 1/10, corresponding to equal changes in 
the densities of the two chains. As Ntot increases fur- 
ther, we see an ellipse (small, red online) moving into 
the rectangle. Here, we might expect the two TV's to in- 
crease together, until the ellipse spans the vertical range 
of the rectangle. From there on, A'2 ceases to change 



Turning next to the study of competition between 
three TASEPs, we see immediately that even more sce- 
narios are possible, from all lengths being equal, to some 
being the same, to all lengths being drastically different. 
Though we have explored quite a few cases, we will only 
present data for the most extreme one (HD): Li = 1000, 
L2 = 100, L3 = 10 with a = 0.7 and /? = 0.3. How the 
three densities vary as we increase Ntot is shown in Fig. 
[T^ We again see the smaller two TASEPs entering into 
a "plateau" regime, with linear profiles and half filling on 
the average. Meanwhile, the longest chain again displays 
a localized shock, as shown in Fig. [131 In more detail, 
the shortest TASEP reaches p = 1/2 first, followed by the 
intermediate chain. The behavior of the longest TASEP 
is much the same as what we observed in the L = 1000 
chain above. It appears that adding a third chain does 
not lead to any qualitatively novel behavior. The overall 
currents, especially if the severe finite size effects associ- 
ated with L = 10 are taken in account, also display few 
surprises. Of course, the phase space for Af = 3 is consid- 
erably larger and new phenomena may very well emerge 
upon closer examination. Further studies are in progress 
and will be reported elsewhere. In the remainder of this 
paper, we will focus on a quantitatively viable, analytic 
picture for an arbitrary number of TASEPs. 



IV. GENERALIZED DOMAIN WALL THEORY 

To understand most of the phenomena we observed, it 
is sufficient to use the simplest approximation [iTj , based 
on self-consistent equations between the feedback depen- 
dence, cteff (Np), and the occupation variables, Ni, for 
the TASEPs. The only serious complication arises in 
the HD case, when aeff reaches (3 and each individual 
TASEP enters an SP-like regime. As presented above, 
a variety of interesting behaviors emerge for which intu- 
itively reasonable, simple arguments paint a good quali- 
tative picture. In this section, we will provide a quanti- 
tative description, which relies on an excellent approx- 
imation, namely, an appropriately generalized domain 
wall theory [13, [3. Proposed about a decade ago for 
the standard single TASEP [ll,[l3, DW theory assumes 
that the configurations are well accounted for by those 
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FIG. 12: (Color online) Average overall density and current 
as a lunction of Ntot for Li — 1000, L2 = 100, and L3 = 10 
with a = 0.7 and /3 = 0.3. 
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with a microscopic interface between two regions, one 
with high density /5+ and another with low density p_. 
The generahzation to a single TASEP constrained by fi- 
nite resources [U [iP'] provided excellent agreement with 
all aspects of simulation data. Referring the reader to 
for details, let us first present a brief summary of 
this approach here, in the context of multiple TASEPs 
competing for the same pool of particles. 

We assume that a configuration of the system of M 
TASEPs can be well approximated by specifying the po- 
sition of the shocks (DWs) on each lattice, ki G [0, L^], 
with £ — 1, . . . , M. Since all chains are subjected to the 
same entry and exit rates of particles, we will further 
assume that the densities before (sites i < k) and af- 
ter (sites i > k) the wall on all chains are identical, i.e., 
p_ = cteff and — 1 — P , respectively. Now, ag// 
depends on the numbers in the pool and so, on the occu- 
pation on the lattices, Ni. Thus, to close the equations, 
we need the relationship between Ni and kf. 



Ne = p^ke + p+{Le - h) (8) 
= {aeff + (3-l)ki + {l-(3)L, (9) 



which leads us to the pool occupation: 
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FIG. 13: (Color online) Average density profile for Ntot = 
600, Li = 1000, 1/2 = 100, I/i = 10 with a = 0.7 and /3 = 0.3. 



and the dependence of cueff on the fc's. In short, from 
aeff = af{Np) = atanh ([TVtot - J2t Ni] N*), we have 



(K) = atanh 



Ntot - (1 - /3) Ei Le - [aeff + - 1) 
N* 



(11) 



r 



Note that Ofe// depends only on the sum 



M 



(12) 



rather than the individual fc's. This simplification will 
be crucial for us to find an exact steady state solution to 
the master equation in our system. Of course, we must 
solve the non-linear, self-consistent equation (fTT|) to de- 
termine cte// (K). This task was performed numerically, 
as in the previous study p^ . Indeed, we have the same 



functional form as before, except that we now encounter 
aeff [K] E^ Li) instead of o-eff (k; L). 

Once aeff (K) is known, the rates for the DW to hop 
to the left {D~) and right (£'"'") can be computed. We 
denote explicitly their if-dependence by a subscript: 



Dk = 



'e//(l - aeff) 
1- P- aeff 

/?(l-/3) 



^ " 1-/3 -a 



eff 



(13) 
(14) 



With these, the master equation for P{{k(} , t), the prob- 
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ability to find DWs at {kg} at time t, is well defined. For 
{kg} lying in the interior of the allowed domain, it is easy 
to write 



M 



dtP{{h},t) 



Y,[DK+,Pm + S,,},t) (15) 

+D+_^P{{ke-Se,},t)] 
-M[D^ + D+] P{{k,},t) 



where S^j is the Kronecker delta. At the boundaries, 
we must impose reflecting boundary conditions on the 
appropriate fc's. Since there are quite a few (3^^ — 1) of 
these conditions, it will be helpful to begin by studying 
the M — 2 case explicitly (4 corners and 4 sides of a 
rectangle). In the next subsection, we will provide some 
details for obtaining the steady-state solution associated 
with this case. We will find that the natural variables 
are {K,Q = ki — k^) rather than (^1,^2)- The insights 
gained here will facilitate the analysis of the arbitrary M 
case, to be presented in the last subsection. 



A. Case with two TASEPs 

Here, we focus on the steady-state solution to Eqn. 
(HH) for M = 2, for < /ci < Li and < /c2 < L2. 
Dropping the t, and setting the left to zero, this equation 
reduces to 

= D-^^^[P{ki + l,k2) + P[ki,k2 + l)] 
+D+_^[Piki - 1, k2) + P(fci, A:2 - 1)] 
-2[D], + D+]P{ki,k2) (16) 



where 



K ^ki + k2 



(17) 



here. The minimum boundary conditions to be imposed 
correspond to one or both DWs being reflected from the 
ends of the lattices. Thus, at the four sides of the rect- 
angle, we write 

= D-^^^[P{ki + l,0) + P{ki,l)] (18) 

+Dl_,Pik, - 1,0) - [D-^+2D+^]P{k,,0) 
= i?-^+jP(0,fc2 + l) + P(l,A;2)] (19) 

+Dl_,P{0, fc2 - 1) - [Dj:^ + 2<]P(0, k2) 
= Dl_^^^_,[P{k,~l,L2) + P{k,,L2-l)] (20) 

+DL.+k,+,Pih + l,L2) 

-\Dt+k,+'^Dl,^+k^P{kuL2) 
= i?++fe,_jP(Li,fc2-l) + P(ii-l,fc2)] (21) 

+^Z,+..+i^(ii.fc2 + l) 



Similarly, the conditions at the four corners are 

= i?r[P(0,l) + P(l,0)]-2i^+P(0,0) (22) 
= i?^,+i,_i[P(Li-l,L2)+P(Pi,i2-l)](23) 

-'^DI^+lP{LuL2) 
= D+^_,P{L, -1,0) + Dl^^,P{L,., 1) (24) 

-\Dt+Dl^]P{L,,Q) 
= Dl^_^P{Q, L2-l)+ Dl^^,P{l, L2) (25) 

-[D+^+Dl^]P{Q,L2) 

However, as discovered in a previous study [l^, there is 
a more subtle boundary condition. For sufficiently low 
iVtot, it is not possible for one or both DWs to reach the 
left boundary. Specifically, if Ntot < (1 — /?) (ii -I- ^2), 
then at least one of the lattices cannot be filled with the 
high density, p+, so that the sum of the DW positions 
must be larger than i^min = ii + L2 — Ntot/ (1 — /?)• 
Another way to understand this limit is that the pool is 
empty {Np = 0) when K reaches i^min- Both ag// and 
vanish and we simply have P = for ki+k2 < ii^min- 
To summarize, we see that P(fci , k2) can be non-zero only 
in a (generally) "cut-rectangular" domain: 

7^ : < fci < Li, 0<k2<L2, ATmin < ki + k2 

(26) 

Taking into account this complication of the boundary 
conditions, and thanks to D being dependent on only 
one variable, these equations can be solved analytically. 
In particular, although the original competing TASEPs 
is a non-equilibrium statistical mechanics problem, the 
DW approximation reduces it to the point that detailed 
balance prevails. Specifically, the product of the rates 
around every elementary loop (through the four points 
(fci, fe) ; (fci + 1, fe) ; (fci + I,k2 + 1); (/ci, + 1)) in 
configuration space is 



^ K+l^ K+l 



(27) 



regardless of the direction taken around the loop. Since 
all closed loops in this space are composed of these ele- 
mentary loops, the Kolmogorov criterion [l^l is satisfied. 
Thus, we have detailed balance 

D^PikiM) = D+_^P{k^-\,k2) = D+-iP{ki,k2-l) 

.(28) 

and the solution can be obtained by recursion. Setting 
P(Li,L2) = ^"' (29) 

where Z is a (normalization) constant, we easily find the 
stationary distribution, for A' < Pi + L2 — 1, 

P{ki,k2) = Z-'^{K) , (30) 

provided {ki,k2) lies in TZ. Here, we define 



^i^)^ n ^ 



(31) 
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and emphasize that P(fci,fc2) varies only via the sum 
ki + k2 and is fiat along lines of constant K. Due to 
the irregular shape of TZ, the normalization factor Z is 
straightforward to find but not simply $ (^); &s in 
the single TASEP case (l9|]. For example, if < fsTmin < 
ii, then 



L1+L2-1 



K) e (Li -K) + L2 + l]^ {K) 

(32) 

where Q is the Heavyside function. 

Since this stationary distribution depends only on the 
sum of the shock positions, it behooves us to change the 
variables from (fci,A;2) to {K,Q), where 

Q = ki-k2 . (33) 

In terms of these variables, Eqn. becomes 

= DJ^^^[P{K+1,Q + 1) + P{K + 1,Q-1)] (34) 
+D+_,[PiK - 1, Q - 1) + P(X - 1, Q + 1)] 
-2[D++Dj,]P{K,Q) 

with similar replacements for the eight boundary condi- 
tions. The detailed balance condition, Eqn. (|28p . now 
reads 



D+P{K,Q) 



Dj,^^PiK + l,Q + l) (35) 
D],^^PiK + l,Q-l) 



and suggests a solution that depends only on K. Assum- 
ing the ansatz 



P{K,Q) = Z-^<i>{K)H(Q,K) 



(36) 



where H (Q, K) is a Heavyside-like function (unity 
in TZ and zero otherwise), we find $ (X + 1) cx 
$ {K) D^/ D~^j^^. Taking some care with the boundary 
conditions, it is easy to verify that the form p6p is indeed 
valid. 

Before turning to a theory for the multi- TASEP case, 
let us compare the predictions (with no adjustable pa- 
rameters) of this approach with the data. First, the den- 
sity profiles can be obtained, following the methods in 
the previous study [l9j : 



Pi{i) = 



L2 

E 

fe2=0 



(1-/?) ^ P(fci,fc2) 



(37) 



fci=0 



E Oieff{K)P{kl,k2] 
fci=i+l 



fei=0 



(1-/3) ^ P(fcl,fc2) 



(38) 



fc2=0 



E "e//(A:)P(fci, ^2) 



/C2=i + 1 



In this expression, it is clear that X^fcs ^('''i' ^2) is just 
the probability for finding the DW at ki regardless of 
the configuration in the other TASEP and similarly for 
J2ki ^'^'^ ihsX these agree well with all observed 

profiles - not only for regimes with little structure, but 
also for complex situations such as the "plateau" region. 
As an illustration, in Fig. [5] we show profiles for the fa- 
miliar Li = 1000, L2 = 100 case with Ntot = 700. It 
is clear that the essentials of our system, such as a lin- 
ear profile in one chain along with a localized shock in 
the other, have been successfully captured in this theory. 
From these profiles, we obtain the overall densities by 
pi — Pi{T')I^i- As shown in Figs. U [3l IH and[5l they 
are in excellent agreement with the simulation data. Not 
surprisingly, the histograms for (A^i, A'2) shown in Fig. [8] 
can also be predicted, being basically P(fci,fc2) via the 
correspondence ([9]). Although the agreement is also quite 
good, we should remark on two shortcomings. First, re- 
lationship ([HI) between iV and k cannot be exact, since 
both N and k are integer- valued. Second, lattices with 
greater than p+L or less than are obviously absent 
from the theory, a limitation due to the DW approxima- 
tion. However, in all the regions we have explored, these 
shortcomings result only in minor disparities. 

Since so many aspects of our system can be understood 
by this approach, let us return to form a better intuitive 
picture for the two "crossover" regions for the L\ = 1000, 
L2 = 100 case discussed at the end of Section llll Dl First, 
note that there are subtle changes in the gradients of Ni, 
as Ntot increases up to the lower crossover (Ntot ~ 400, in 
Fig. [11]). By contrast, dNp/dNtot remains relatively con- 
stant. More significantly, the longer chain "loses" while 
the shorter one "wins." This behavior can be traced to 
the DWs being mostly bound to the exit (right edge), but 
making longer and longer excursions into the lattice as 
particles in the system become more abundant. The ex- 
ponential tails of this excursion are essentially identical, 
provided they do not intrude significantly into the entry 
(left edge). Indeed, if we measure the ratio of the profiles 
pi (900 -I- i) I p2 (i), it is essentially unity for Ntot up to 
~ 400. Nonetheless, the overall density of the shorter lat- 
tice is affected more by a similar portion of an enhanced 
profile (over p_), and so, will increase faster. We should 
also remark that the density to the left of the shock (low 
density region) in either lattice has yet to reach the final 
value (i.e., p_ < /3). 

After we enter the first crossover region {Ntot ~ 400), 
the DW wanders further from the exit in each lattice. 
It eventually reaches the left side of the smaller TASEP 
and enters the SP. To understand how this transition oc- 
curs, let us look at how the density of the ordinary, un- 
constrained TASEP of infinite length changes as a func- 
tion of a with P < 1/2, as we cross the first order line 
a — f3. Starting from a = 0, the density rises linearly un- 
til a = (3. At this point, the density jumps to 1/2. It then 
jumps again to a value of 1 — /3 for a> (3. For a system 
with a finite length, these jumps are no longer sharp, but 
smeared out near the value of a = (3, reflecting a rapid 



12 



increase in the density (and number of particles) in this 
region. We now use this information, noting that (i) Nfot 
controls the effective a, and (ii) ae// = P spans a whole 
region of Ntot, and (iii) that the shorter TASEP behaves 
essentially like an unconstrained one, since it can draw 
from both the pool and the longer chain. The first rapid 
increase, as aeff approaches f3 from below, is responsi- 
ble for the first peak of dN2/dNtot just above Ntot = 400, 
shown in Fig. [TTJ Clearly, the gradient then decays to 
zero as N2 reaches the characteristic plateau. The next 
peak reflects the end of the ag// = (3 region: it is related 
to the sharp increase of the density on the other side of 
the first order line. In terms of the domain wall picture, 
the first peak reflects the fact that the exponential tail in 
the density profile reaches the entrance and changes from 
an exponential decay into a linear one. In other words, 
the probability of finding the shock near the entrance in- 
creases until it is flat across the whole system. At the 
second crossover, the reverse happens: the system is in a 
high-density phase, with a small exponential tail at the 
entrance, indicating that the shock is now predominantly 
found there. 

It is interesting to note that, in the first crossover 
regime, the pool "loses" steadily at sharing the increases 
in Ntot- Of course, over the plateau regime, only the 
longer chain gains from the changes in Ntot ■ The picture 
for the second crossover regime is essentially the same, 
except occurring in "reverse order." Needless to say, the 
pool is the ultimate "winner" in this competition, ab- 
sorbing all increases in Ntot beyond this regime. 

B. Multiple TASEPs 

The insights gained in the detailed analysis for two 
TASEPs greatly facilitate investigation of the general 
multi- TASEP case here. In particular, since the entry 
rates into all M chains are the same, cteff, and depend 
only on K = kg, the transition probabilities in the 
master equation again satisfy the Kolmogorov crite- 
rion. A change of coordinates similar to the two TASEP 
case can be performed, with i^T = ^ fc^ as the special 
variable. The boundary conditions are also straightfor- 
ward though considerably more tedious, since the M- 
dimensional generalization of TZ is more complicated in 
general. Nevertheless, thanks to detailed balance, it is 
simple to solve for the stationary distribution associated 
with this (effectively one-dimensional) master equation. 
The answer will again be the form (j36p : 

P{{ke}) = Z-^<S>{K)H{{k,}) (39) 

where the Heavyside-like function is now defined for the 
M-dimensional TZ. Needless to say, the "cuts" of constant 
K across TZ are geometric objects of dimension Af — 1, 
the general shapes of which are far more complicated 
than the lines in the M = 2 case above. As a result, 
the normalization factor is given by a much more com- 
plex expression than (|32p . since the coefficients involve 



polynomials (in K) up to order M — 1. Nevertheless, 
these computations are very simple for modern comput- 
ers and P({fcf}) can be easily accessed numerically, for 
a reasonable range of Ms. For example, we computed 
P{ki,k2,k3) for the three TASEP case discussed in the 
previous section {Li = IOL2 ~ IOO-L3 = 1000, with 
a — 0.7 = 1 ~ (3). From P, we obtained the average 
profiles and overall densities, as a function of Ntot- As 
shown in Fig. I12i there is excellent agreement between 
this theory and the data, over the many regimes encoun- 
tered while Ntot increases. Similarly, we find good agree- 
ment for the average density profiles, as illustrated in 
Fig. [13] for the case with Ntot — 600. While two profiles 
are linear and one contains a localized shock, all proper- 
ties are well predicted by this approach. Of course, the 
crossover regimes are richer, as each of the shorter chains 
produces peaks in the gradients similar to those shown 
in Fig. 111! With no qualitatively new phenomena, these 
and other cases, as well as further details, will be dis- 
cussed in a later publication [25| . Our conclusion is that 
the generalized DW theory is remarkably successful at 
capturing the essence of our problem. Only the study of 
very sensitive quantities reveals a poorer match between 
this approach and simulation data. 

To end this section, let us note that there is a ther- 
mal analog for P({ki}), namely, the canonical ensemble. 
Here, K plays the role of the total energy, E, while the to- 
tal number of points in a sheet of fixed K corresponds to 
the microcanonical partition function, Q (i?). The point 
{ke} — {Li} (highest K allowed) would be the "ground 
state," while the precise connection between E and K 
is given hy E = — fesrin^ {K). Of course, we chose to 
label our normalization constant Z to carry this anal- 
ogy to its logical end. Meanwhile, Ntot seems to play 
the role of temperature, with the average of —K being a 
monotonically increasing function of Ntot - Exploring this 
correspondence will be both interesting and imperative, 
especially if we hope to make progress toward our goal, 
namely, a cell with a few thousand different genes, each 
appearing in hundreds of copies. 



V. SUMMARY AND OUTLOOK 

In this paper, we explored how competition between 
TASEPs affects the density profile, the overall density, 
and the current for each chain. A feedback mechanism 
introduced previously [l3| was implemented. We used 
Monte Carlo simulations to explore the properties of the 
overall densities, currents, and profiles of the TASEPs, 
for a variety of parameters (lattice length, Ntot, ct, and 
(3). The competition produced several novel features that 
are absent from the case of a single TASEP constrained 
by finite resources [13, Ull- There, the feedback serves 
to localize a domain wall, when the control parameters 
are set in favor of its appearance in the lattice. Here, 
the presence of other TASEPs adds an extra dimension 
to the feedback, leading to the delocalization of the DW 
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to the extent allowed by the lengths of the other chains. 
Thus, when a long chain competes with a short one, its 
DW wandering is limited by the length of the latter. By 
contrast, a DW in a short chain is free to roam over 
the entire lattice, so that the average profile displayed is 
strictly linear and the average overall density is just 1/2. 
This picture can be readily generalized to three or more 
chains and is confirmed in limited simulation studies with 
three competing TASEPs. 

For the single TASEP with feedback, the standard 
domain wall theory was appropriately generalized and 
proved to be extremely successful 19]. Extending this 
theory to an arbitrary number of TASEPs is straight- 
forward and a master equation for P {{ki} , t), where ki 
denotes the position of the DW in the chain, is easily 
formulated. Fortunately, we are able to find the steady- 
state solution in a system where the entry rate onto all 
chains arc the same. The details were presented for the 
two TASEP case and steps for extending it to arbitrary 
M were provided. A remarkable features is the existence 
of an intimate mapping from our steady state to a canon- 
ical equilibrium ensemble. From this stationary distribu- 
tion, all density profiles and currents can be predicted, 
with no adjustable parameters. The results agree well 
with all data and give us much insight into a variety of 
phenomena discovered first in simulations. 

The most intriguing questions beyond our study here 
focus on time-dependent phenomena. Even for the stan- 
dard TASEP, there is a wealth of interesting dynamics 
0, How are these affected by constraints of finite 
resources and competition? For example, one of the sim- 
plest quantities displaying remarkable behavior in the 
open TASEP is the power spectrum associated with the 
total occupation, N (t) Preliminary data for a single 
chain coupled to finite resources reveal a host of novel 
phenomena [26| . With many chains competing for one 
pool of particles, we can study many other quantities, 
such as correlations between the various TASEPs. Hope- 
fully, these pursuits will reveal other exciting secrets in 
this system and provide us with a deeper insight. 

Being motivated by protein synthesis in cells, there are 



many extensions we can explore. Here we list some exam- 
ples. The initiation rates for various genes are far from 
being identical. Thus, we should introduce a full set of 
{a^}s, to model highly and rarely expressed genes. Natu- 
rally, these will induce complex {fc^}-dependences in the 
DW hopping rates: ({fc^}). One serious consequence 
is that the transition rates in the master equation (jl6|) 
will typically violate the Kolmogorov condition, so that 
the steady-state solution will be truly non-equilibrium in 
character. Nontrivial steady-state probability currents 
necessarily follow [23l |. and their implications surely de- 
serve further pursuit. Other obvious extensions include 
important aspects of protein synthesis which have been 
considered in earlier models, such as having particles with 
finite extent to model the fact that ribosomes are rela- 
tively large molecules "covering" many codons [3], and 
the inclusion of inhomogeneous hopping rates along the 
lattice to model the inhomogeneous sequence of codons 
and the wide range of concentrations of their associated 
aa-tRNAs [13] . Another aspect is the ribosome recycling 
enhancement considered by T. Chou [l3| and the role of 
diffusion of (the subunits of) ribosomes in a competitive 
environment. Along these lines, to model the workings of 
a cell better, we should consider a system with some reg- 
ulation on Ntot, as opposed to being just a preassigned, 
fixed number. Finally, an even more ambitious goal is to 
include not only the competition for ribosomes, but also 
for the many varieties of aa-tRNA molecules. Clearly, 
much work remains to be done in order to arrive at a 
realistic model of, and to better understand, protein syn- 
thesis in a cell. 
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